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SUMMARY 

The modal element method for acoustic scattering can be simplified when the scattering body is rigid. In 
this simplified method, called the modal ring method, the scattering body is represented by a ring of triangular 
finite elements forming the outer surface. The acoustic pressure is calculated at the element nodes. The pressure 
in the infinite computational region surrounding the body is represented analytically by an eigenfunction expan- 
sion. The two solution forms are coupled by the continuity of pressure and velocity on the body surface. 

The modal ring method effectively reduces the two-dimensional scattering problem to a one-dimensional 
problem capable of handling very high frequency scattering. In contrast to the boundary element method or the 
method of moments, which perform a similar reduction in problem dimension, the modal line method has the 
added advantage of having a highly banded solution matrix requiring considerably less computer storage. The 
method shows excellent agreement with analytic results for scattering from rigid circular cylinders over a wide 
frequency range (1 = ka 2 100) in the near and far fields. 


INTRODUCTION 

The modal element method, which couples finite elements and eigenfunction expansions, has been em- 
ployed in electromagnetic and acoustic scattering and duct transmission problems. The primary reasons for 
employing this technique are (1) to describe accurately the radiation boundary condition at the computational 
boundary and (2) to reduce the size of the numerical grid. This hybrid steady state method has been given 
various titles, such as the unimoment method, the transfinite element method, and the modal element method. 

In electromagnetics, Chang and Mei (1976) and Lee and Cendes (1987) applied the method to scattering from 
dielectric cylinders while Baumeister (1991) applied the method to electromagnetic propagation in ducts. In 
acoustics, Astley and Eversman (1981) employed the method in duct propagation problems while Baumeister 
and Kreider (1993) have applied the method to acoustic scattering problems. 

The purpose of this paper is to develop and study an extension of the modal element method for scattering 
from rigid bodies. The goal is to minimize the domain in which finite elements are employed. This approach is 
called the modal ring method and was briefly introduced by Baumeister and Kreider (1993). 

The modal ring method can effectively reduce a two-dimensional scattering problem to a one-dimensional 
problem by employing a ring of nodes along the surface of the rigid scattering body. Herein, a variety of grid 
configurations and parameters are explored. The method is validated by applying it to the simple case of scat- 
tering from a hard circular cylinder, for which analytic solutions are easily obtained. 
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NOMENCLATURE 


A n) + modal amplitude of wave moving in +r direction away from origin 
a dimensionless circular cylinder radius 

e r relative error 

H m ^ Hankel function of the first kind 

i fT 

k wave number 

m mode number 

M CO ef number of modal coefficients used in eigenfunction expansion 

p dimensionless perturbation acoustic pressure 

r dimensionless radial coordinate 

e dimensionless complex acoustic permittivity 

0 angle between position vector and x axis 

p dimensionless complex acoustic permeability 

to dimensionless frequency 


METHOD OF ANALYSIS 

This study is concerned with computing the acoustic scattering by a two-dimensional rigid body of an 
impinging plane wave traveling in the +x direction. The spatial domain is divided into two subdomains, the 
finite element domain, which contains the body, and the homogeneous domain, which surrounds the body and 
extends to infinity. Linear triangular elements are used in the finite element 

domain to calculate the pressure at the nodes. In the homogeneous domain, an eigenfunction expansion repre- 
sents the acoustic pressure. The two solution forms are coupled by imposing continuity on the pressure and 
velocity at the interface between the two subdomains. This coupling results in a single matrix equation in which 
the eigenfunction coefficients and the pressures at the finite element nodes are calculated simultaneously, yield- 
ing a global representation of the acoustic field. 


GEOMETRICAL MODEL 

For penetrable bodies, the internal pressure field must be known to determine the scattered field. The 
modal element method (Baumeister and Kreider, 1993) uses a finite element grid inside the body to calculate 
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the internal pressure Held. For rigid bodies, studied here, the field does not penetrate the body, so the grid can 
be greatly reduced, yielding a variation of the method called the modal ring method. 

In the modal ring method, the scattering body is represented by a thin ring of finite elements following the 
body’s contour. Figure 1 shows three possible configurations. The grid in figure 1(a) consists of a single row of 
boundary nodes with a central node. The grid in figure 1(b) has a ring of elements on the body’s surface and a 
central node, so the body is filled with elements while the grid in figure 1(c) has the same ring of elements but 
excludes the central region, where the field vanishes. The rigid body is simulated numerically, through an impe- 
dance mismatch induced by setting e = 1 - 10 19 i and p = 1 at each internal finite element node. This feature 
allows greater flexibility in the numerical implementation of the method, because penetrable or coated rigid 
bodies may be studied with only slight modifications to the computer code, mainly in grid generation. For pene- 
trable bodies, the internal grid is more extensive, while for coated bodies, several rings in the coating region 
may be needed. 


GOVERNING EQUATIONS 


Acoustic propagation in two-dimensional space can be modeled by the continuity, momentum, and state 
linearized gas dynamic equations in the absence of flow. For harmonic pressure propagation in an inhomoge- 
neous material, the following dimensionless wave equation applies: 


a 

f 1 d P ) 

a 

+ 

f 1 dp ' | 

dx 

i e8 *j 

dy 

l* 


to 2 | ip = 0 


( 1 ) 


The harmonic time dependence e' iwt has been factored out. In analogy with electromagnetic scattering, e is the 
acoustic “permittivity," p is the acoustic “permeability" and to is the dimensionless frequency. The wave num- 
ber is 

( 2 ) 

k = \Zto 2 pe 


At the interface between the finite element region and the analytic region, continuity is imposed on the 
pressure and velocity. The radiation boundary condition at infinity is automatically satisfied by the eigenfunction 
expansion introduced in the next section. 


ANALYTIC SOLUTION 

In the homogeneous domain, an exact eigenfunction expansion can be derived from equation (1) by 
separation of variables. An incident plane wave travelling in the +x direction that strikes a symmetric two- 
dimensional scatterer generates a pressure field approximated by 

A i coe f\ 

p = pi +p s . e ikx + g A^H^\kr)cos(md) 


p* is the incident plane wave; p s is the scattered wave. The modal coefficients A m + are unknowns to be deter- 
mined. Formulas to estimate the number of modes needed for convergence can be found in Baumeister and 
Kreider (1993). 
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FINITE ELEMENT SOLUTION 


The finite element domain is divided into triangular elements with unknown acoustic pressure at the nodes. 
It is assumed that all material properties are constant in each element Details on die numerical setup and solu- 
tion can be found in Baumeister and Kreider (1993). 


RESULTS AND COMPARISONS 

The modal ring method may be implemented using various grid systems and parameter values for ring 
thickness and the number of nodes. In order to study the influence of these parameters on numerical results, the 
following simple standard scattering problem is considered: a unit plane wave, incident from the left with a fre- 
quency ranging from k = 0.1 to k = 100, strikes a rigid circular cylinder of dimensionless radius r = 1 oriented 
with its axis normal to the propagation direction. In the homogeneous region surrounding the scatterer, e = 1 
and p = 1, while within the scatterer, e = l-10 19 i and ji = 1. The number of modes M^ used in the eigen- 
function expansion in equation (3) is determined by formulas in Baumeister and Kreider (1993, eqs. (26) 
and (27)). All of the plots presented here are generated using equation (3) with the calculated A m + values. 


Number of Nodes 

As the frequency of the incident wave increases, more finite element nodes are required to resolve the 
increasing oscillations in both the pressure and intensity fields. The oscillations in this case are represented by 
the cosine terms in equation (3), so the number of boundary nodes depends on the number of eigenfunction 
modes used. Experience has shown that setting the number of boundary nodes to ^M^f suffices; this condi- 
tion is used in the examples shown below. 


Grid Systems 

Grid 1(a), with its single ring of boundary nodes, is attractive because it uses the least number of nodes for 
a given problem. However, it directly links the boundary of the body with the center node. This inflexibility 
yields poor numerical results — relative errors up to 30 percent for low frequencies. This grid should not be used 
for any scattering problem. 

In contrast, grids 1(b) and (c) consist of a ring of elements following the body’s surface. Grid 1(b) con- 
tains a central node, while grid 1(c) excludes the central region. When using grid 1(c), the pressure gradient 
along the inner element boundaries is set to zero. Numerical results for these grids are virtually identical, indi- 
cating that the presence of the central node is not critical. Figure 2 shows the pressure amplitude at radius 
r = 3.18 with frequency ka = 25 using grid 1(b). Here, M coef = 40 modes are required for the scattered field in 
equation (3), so 480 boundary nodes are used. Figure 3 shows the pressure amplitude at radius r = 3.18 with 
frequency ka = 50 using grid 1(c) with 804 nodes (M coef = 67). In both figures, the numerical results (squares) 
are in excellent agreement with the analytic solutions (solid lines). 

Both grids were tested with frequencies up to ka = 100 with accuracy comparable to that seen in figures 2 
and 3. The plots are not included here because the high degree of oscillation in the graphs makes visual inspec- 
tion difficult. 
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Ring Thickness 


The thickness of the ring of elements surrounding the body is very significant. For a cylinder of radius 1, 
element thicknesses ranging from 0.05 to 0.001 were considered. Figure 4 shows results for incident frequency 
k = 25. Here, M^f = 40, so 480 boundary nodes are used. Grid 1(b) is used. For element thickness 0.05 
(fig. 4(a)), results are poor. For element thickness 0.01 (fig. 4(b)), results are better, but there is still some 
deviation from the exact solution. Reducing the element thickness to 0.001 (fig. 4(c)) yields excellent results. 


CONCLUDING REMARKS 

The modal ring method for acoustic scattering from a two-dimensional rigid body is presented. The acous- 
tic pressure field is represented by a finite element solution on a thin ring of elements on the surface of the 
scattering body, and by an eigenfunction expansion outside the body. The two representations are coupled by 
the continuity of pressure and velocity across the interface between the two subdomains, and are calculated 
simultaneously from a single matrix equation. This matrix is highly banded because the usual two-dimensional 
finite element grid is replaced by an essentially one-dimensional grid. This reduction in computational resources 
is the main benefit of the modal ring method, but it may be used only when the scattering body is rigid. How- 
ever, the method is applicable to problems involving very high or low frequency scattering from rigid bodies of 
complicated shape. In the validation cases presented, numerical results are in excellent agreement with the cor- 
responding exact solutions. 

The grid system most suitable for use with the modal ring method contains a very thin ring around the 
surface of the scattering body, and may or may not include a node at the body’s center. If the central node is 
not included, the pressure gradient along the inner element boundaries is set to zero. 
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Exact analysts ABS (p) 

□ Numerical solutions ABS (p) 



Figure 2.— Polar Plot of the total acoustic pressure around a rigid 
cylinder modeled by a single line of symmetrical elements on 
boundary with central node for configuration figure 1 (b) (ka = 25,0, 
a = 1, number of nodes = 961 , number of elements = 1440, 

Mcoef = 40, At = 0.001 , r = 3.1 8, kr = 79.58). 


Exact analysis ABS (p) 

□ Numerical solutions ABS (p) 



Figure 3. — Polar Plot of the total acoustic pressure around a rigid 
cylinder modeled by a single line of symmetrical elements on 
boundary for configuration figure 1 (c) (ka = 50.0, a = 1, number 
of nodes = 1 608, number of elements = 1 608, M COO f = 67, 

At = 0.001, r = 3.18, kr = 159.1 6). 
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(b) At = 0.01, (c) At = 0.001. 

Figure 4. — Sensitivity of the total acoustic pressure around a solid cylinder to line element thickness At (ka = 25.0, a = 1 , r = 1 .5, 
number of nodes = 961 , number of elements = 1440, M coef = 40, configuration figure 1 (b». 
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